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Abstract. The equation which describes a particle diffusing in a logarithmic 
potential arises in diverse physical problems such as momentum diffusion of atoms in 
optical traps, condensation processes, and denaturation of DNA molecules. A detailed 
study of the approach of such systems to equilibrium via a scaling analysis is carried 
out, revealing three surprising features: (i) the solution is given by two distinct scaling 
forms, corresponding to a diffusive {x ^ y/t) and a subdiffusive {x ^ length 
scales, respectively; (ii) the scaling exponents and scaling functions corresponding to 
both regimes are selected by the initial condition; and (iii) this dependence on the 
initial condition manifests a "phase transition" from a regime in which the scaling 
solution depends on the initial condition to a regime in which it is independent of 
it. The selection mechanism which is found has many similarities to the marginal 
stability mechanism which has been widely studied in the context of fronts propagating 
into unstable states. The general scaling forms are presented and their practical and 
theoretical applications are discussed. 



PACS numbers: 05.40.-a,05.10.Gg 
1. Introduction 

A large variety of physical problems are governed by the simple diffusion equation which 
describes a Brownian particle in a logarithmic potential. Such problems range from the 
momentum spreading of cold atoms in optical traps to the dynamics of "bubbles" 
in denaturing DNA molecules [6l-lll[. and from the relaxation of a single particle in a 



fluid with long-range interactions |12|-|14| to models describing brief awakenings in the 
course of a night's sleep 15|. In these examples, as well as in others which will be 
described below, one is interested in the distribution of a fluctuating quantity x: for 
optically trapped cold atoms x stands for the momentum of the atom, in the problem of 
DNA denaturation x is the length of a denatured unbound loop in the double stranded 
molecule, and when modeling sleep dynamics, x represents the wakefulness level of 
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a sleeping individual. In the problems we consider, the temporal evolution of the 
distribution P{x,t) can be approximated when x is large enough by the equation 

dPi^.t)_a,,,.^rm^^ (1) 



dt dx 



X 



dx"^ 

The dimensionless parameter 6, which plays a central role in the solution of the equation, 
has a different physical meaning in each problem. In physical systems of the type 
which we consider, Eq. ([T]) has corrections at small values of and, in particular, the 
divergence at small x is not present. 

To have a concrete physical picture in mind, we will mostly concentrate on a 
Brownian particle diffusing under the influence of a one- dimensional external potential 
which is logarithmic for large x: 

^^^^Mog(W). (2) 

KbT 

where T is the temperature, fee is Boltzmann's constant and h measures the strength of 
the potential. The corresponding Fokker-Planck equation has the form of a continuity 
equation for the probability distribution 

^ = -^. with = (3) 

where J{x, t) denotes the probability current. To simplify notation, here and throughout 
the paper time is measured in units in which the diffusion coefficient is equal to 1, and 
the potential is measured in units of temperature, i.e., ksT = 1. For large values of |x| 
this Fokker-Planck equation reduces to ([1]). For some applications it is natural to restrict 
the variable x to be non-negative, in which case the equation should be supplemented 
by a boundary condition at x = 0. We therefore discuss both the case of restricted x 
with different boundary conditions, and the case where x is unbounded (which requires 
no additional boundary conditions). 

Our goal is to study the long-time behavior of the solutions to this diffusion problem. 
In a recent paper [l6| we have reported that the solutions of equation ([3]) relax to 
equilibrium via a universal scaling form which depends on the potential only through 
its asymptotic form ([2]), and also depends on the initial condition. In the present paper 
we elaborate on the analysis of 16|, present in more detail the derivation of this result, 
and discuss its implications and applications. 

The scaling form which we find exhibits several features which are not typically 
found in scaling solutions, (i) For any large finite time, the overall scaling form is 
comprised of two scaling functions, one corresponding to small values of x {x <ti 
and another corresponding to large values of x (x ~ y/t). Together, these two scaling 
functions give the distribution P{x, t) for all values of x, including the microscopic scale 
where the potential deviates from a logarithm. The small-x details of the potential 
V{x) enter only through the steady-state distribution which multiplies this scaling 
form, (ii) The equation admits not only one, but a family of such overall scaling forms, 
each characterized by different scaling exponents. The scaling form which describes 
the observed relaxation to equilibrium depends on the initial condition via a selection 
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mechanism akin to the marginal stabihty selection mechanism encountered in, e.g., 
fronts propagating into unstable states, (iii) As in problems of propagating fronts, by 
continuously changing the tail of the initial condition, the selected scaling exponent 
exhibits a "phase transition" from a smoothly varying to a fixed value. 

Features (ii) and (iii) of the scaling solution provide an interesting connection 
between this diffusion problem and the well-known problem of fronts propagating into 
unstable states. Many systems of the latter type admit a family of traveling-wave 
solutions with different propagation velocities, and the mechanism by which the eventual 
velocity and waveform are selected has been widely studied As described below, 
the selection mechanism which we find for the diffusion problem is similar in many of its 
details to that corresponding to propagating fronts. The two problems differ however 
in some basic aspects: unlike the homogeneous nonlinear propagating front problems, 
equation (jS]) is linear yet inhomogeneous in space. Although the inhomogeneity of our 
problem restricts the utility of mathematical methods used to analyze the selection of 
front velocities, most notably Fourier analysis, its linearity makes it exactly solvable and 
facilitates the demonstration of the selection mechanism. The similarities between the 
problems, suggests that a common mathematical description of their solutions might 
exist. 

As mentioned above, at late times the entire distribution is given by a scaling form 
(feature (i) above). At late times t ^ 1, two different length scales emerge: a large-x 
length scale of a; ~ t^^^ and and a small-x length scale of x ~ t"^ with a 6-dependent 
exponent 7 < 1/2. Hereafter we refer to these two length scales as "the large-x" and "the 
small-x" regions, respectively. The exponent 7 depends on the boundary conditions, 
and in particular, for a reflecting boundary condition at the origin, 7 = 1/(6 + 1). 
The solution in each of the length scales is given by a different scaling function, with 
a smooth interpolation between the two functions. Moreover, to leading order in t, 
these two scaling functions yield the solution at any point x (see Fig. [1]). Both scaling 
functions are selected by the initial condition: the large-x scaling function determines 
the one in the small- a; region. In the language of traveling waves, this corresponds to 
a system with two fronts propagating with different velocities, whereby the selected 
scaling solution of the "faster front" dictates that of the "slower" one. 

For a wide class of initial conditions, which include compactly supported (or 
"steep") initial distributions, the large- a: scahng solution of Eqs. ([2])-([3]) has recently 



been found in p^, and [20[. There, the dependence of the solution on the initial 
condition and the behavior at small x have not been considered. In fact, in many 
physical circumstances the initial distribution is not steep. In other situations, the 
small-x behavior rather than the large-x one determines the physical observables of 
interest. In these two respects, beyond the relevance of our work to the general theory of 
selection problems and scaling solutions, it also presents a comprehensive analysis of the 
scaling solution of (I2])-(|3]) which provides useful results for many concrete systems. To 
demonstrate the applicability of our results, we consider at the end of this paper three 
physical examples: (a) cold atoms in an optical lattice undergoing a rapid "quench" 
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from one steady state to another. Here we discuss initial conditions with a fat tail, (b) 
Nonequilibrium driven models exhibiting real-space condensation. Here we show that 
current correlations in these systems may be evaluated by considering initial conditions 
with specific algebraic decay at the tails, (c) The dynamics of loops in DNA molecules 
undergoing denaturation. Here the effect of an absorbing boundary condition is probed. 

The paper is organized as follows. In Sec. |2] we present the scaling solution and 
its selection mechanism, and heuristically derive its form. The results of this section 
are backed up by an exact solution of the Fokker-Planck equation ([T]) which appears in 



Appendix A In Sec. |3], the scaling solution is discussed in the broader contexts of the 
general theory of scaling solutions (Sec. 13. ip . of selection in problems of propagating 
fronts (Sec. 13. 2p . and of the results of previous work on Eq. ([1]) (Sec. 13. 3p . The 
discussion in Sections |2] and |3] is focused on systems whose boundary conditions conserve 
probability. In Sec. |4] we generalize the results of the previous sections to the case in 
which probability is not conserved at the boundary. In particular, we show that the 
large- 1 a; I scaling form is not affected by the boundary conditions. Applications of our 
results to concrete physical systems are discussed in Sec. [5], in which we also present a 
general review of some of the problems which are described by Eqs. (I2])-(I3]), for which 
our results may be relevant. Finally, Sec. [6] contains a summary of our results and some 
concluding remarks. 

2. The scaling solution and its universal character 

In this section we present the scaling solution of Eq. (I2])-(|3]). We begin in Sec. 12. II with a 
general discussion of the problem of diffusion in a logarithmic potential, and present its 
scaling solution. In the following subsections this results is derived heuristically, while 



the exact derivation of this result, which is somewhat technical, is found in [Appendix A 
First, in Sec. 12.31 we demonstrate that in the large-x regime of x ~ -\/t, Eq. ([1]) admits a 
one-parameter family of scaling solutions. In Sec. 12.41 we present the selection criterion 
which explains how the initial conditions determine which member of this family is 
eventually observed. In Sec. 12.51 we derive the scaling form for the small- a; regime. The 
derivation of Sec. 12.31 rests on the assumption that when x y/t ^ 1, Eq. (|3]) can be 
well approximated by Eq. ([1]). In Sec. I2.6l we justify this assumption by showing that our 
scaling solution is universal, i.e., it depends only on the large-|a;| tails of the potential 

2.1. General discussion of the problem 

We begin by writing the Fokker-Planck equation (ISJ-Q in more concrete terms. We 
consider a potential has of form 

V{x) = b\og{\x\) + U{x), (4) 
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where the correction U{x) is neghgible for large |x| and it ensures that V{x) does not 
diverge at the origin. For concreteness, we assume that for large x 

U{x > 1) = 0(|s|-'^) with (7 > 0. (5) 

Regularizing the potential at small x is needed since for 6 > 1, the case on which we 
focus below, a logarithmic divergence of the potential at the origin makes x = an 
absorbing state and any normalized initial condition tends to a (5-function distribution 



around it [21|, |22|. This suggests that in systems with 6 > 1, physical corrections to the 
logarithmic potential near the origin cannot be neglected when analyzing the long-time 
behavior. 

With this notation, the Fokker-Planck equation (j3]) is 



dP{x,t) d r6,^ , , ^, ,] d^P{x,t) 



dt dx 



X 



l + h{x))P{x,t) +—^^ (6) 



where 



xU'(x) 

h{x) ^ = 0{\xr). (7) 

We begin in this section by considering only the boundary condition at x = where the 
probability flux at the origin is zero, i.e., J(0) = 0. This corresponds to diffusion on the 
entire real line with an even initial condition, or diffusion on the positive half line with 
a reflecting boundary condition at the origin. The general case, and the effect of other 
boundary conditions will be examined in Sec. HI 

The stationary solution of the diffusion equation ([3]) in this potential has the form 
of a Boltzmann distribution 

P*{x) = |e-^(^) ~ ia:-^ (8) 

where Z is a normalization constant given by Z = / e-^(^)cix. For b > 1, Z is finite 
and the system tends towards this unique equilibrium distribution regardless of the 
initial condition (we assume that the potential does not contain infinite energy barriers 
and the system is ergodic). However, for b < 1, the equilibrium distribution cannot 
be normalized. In this case, any normalized initial condition tends to zero. Thus, 
potentials with logarithmic tails are a marginal case for the diffusion equation. Any 
potential which increases at large x faster than a logarithm "traps" the particle and 
the probability distribution reaches a steady state at long times. On the other hand, 
potentials which increase with x slower than logarithmically are non-trapping and the 
probability distribution eventually spreads out to infinity. In the marginal case where 
the potential is logarithmic at large x, the particle is trapped at low temperatures and 
becomes delocalized at high temperatures, as the dimensionless parameter b changes 
from 6 > 1 to 6 < 1. 

The aim of this paper is to describe how P{x, t) relaxes towards the eventual 
Boltzmann distribution. We therefore concentrate on the normalizable case 6 > 1. 



Recently, 18N20| have shown that this relaxation is given by a useful and compact 



scaling form. The scaling form which was found, however, describes the solution of 
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Eqs. ([S])^© only for a specific (albeit large) class of initial conditions, and it represents 
correctly only the large- |a;| behavior of the actual scaling solution. The main result of 
the present work is the surprising fact that the long-time scaling form of the solution 
depends on the initial condition in a non-trivial fashion. Moreover, the entire solution 
can be described by a scaling form, where the non-universal features are contained in 
P*{x). We now present these results, and then derive them. 

As we are interested in the relaxation towards the equilibrium distribution, it is 
convenient to study the deviation of P(x, t) from P*{x). To this end we define a function 
G(x, t) via 

P(x,t) = P*{x)[l + G(x,t)], or equivalently G(x,t) = ^^'^'p.^^j'*^'^^ • (9) 

We seek scaling solutions of P{x,t) — P*{x), or equivalently of G{x,t), rather than 
of P{x,t). Note that since the Fokker-Planck equation is linear and it is satisfied by 
P*{x), the distribution P{x,t) and the deviation from equilibrium P{x,t) — P*{x) = 
P*{x)G{x,t) satisfy the same equation. However, while J P{x,t)dx = 1, here 
J P*{x)G{x,t)dx = 0. Looking for a scaling form for solutions with zero normalization 



is our primary extension of the calculations of |18|-|20| which enables us to find all scaling 
solutions to the problem. 

2.2. The scaling solution 



The Fokker-Planck equation ([T]) can be solved exactly by standard methods [23j. By a 
transformation of variables, it can be mapped to a Schrodinger equation in imaginary 
time which describes a quantum particle moving in an inverse square potential. Analysis 
of this quantum problem yields the exact solution of the equation for arbitrary initial 
conditions represented as a series of Bessel functions. Asymptotic analysis of these 
Bessel functions allows one to identify the scaling form which characterizes the approach 
to equilibrium at late times. Although this calculation is straightforward, it is rather 



technical and lengthy. We therefore delay its presentation to Appendix A Here, we 
present the results of this calculation, and in the rest of this section we heuristically 
motivate these results. 



According to the exact calculation of Appendix A the long-time behavior of the 



solution of the Fokker-Planck equation ([3]) with any normalizable initial condition and 
for a reflecting boundary condition at the origin, is given, to leading order in t, by 

{gfi{-^rMn))t~^ for kl < xi{t) 
(10) 
/^(^)r^ for \x\ > x,{t) 

where xi{t) can be chosen to have any value which satisfies 

tV(^+l)«Xi(t)«tl/2, (11) 

and the scaling functions are given by 

— 4(&+l) I ^b+l 
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Figure 1. A schematic representation of the sohition P{x, t) (Eq. ((TO)) ) at a given late 
time i 3> 1 (not drawn to scale). The red double line represents the small-x scaling 
form gj3{x/t^^''^^'^^)t~^ while the blue solid line represents the large-x scaling form 
fp{x/t^^'^)t~^ . The interval on which the two solutions overlap (Eq. (fTT|)) is labeled 
xi{t). 



/,(«) = «^+^F,(l±^;^;-^j. (13) 

The values of the scahng exponents (3 and 6 and of the constant C will be discussed 
shortly. The constant r depends on the domain on which the diffusion is defined: 
r = 2 for diffusion on the positive half-line (with a wall at the origin), while r = 1 for 
symmetric diffusion on the entire real line. Here, iFi is the confluent hyper geometric 



function, whose known properties yield the asymptotic form 24 

for M < 1 

for m > 1, /3 < 1 , (14) 
M^+^e""^ forM>l,/3 = l 

where D = ^'^'^''njb+m] 
wneie r(/3-l) 

The solution (ITOl) is presented schematically in Fig. [H It is made up of two different 
scaling forms with different dynamical exponents: x ~ t^/C'+i) and x ~ t^^^. At small 
values of \x\, the solution is fiat up to x ~ t^^^''~^^\ with a value that approaches zero as 
t~^. At large values of \x\, the solution exhibits a peak at x ~ y/t, whose height shrinks 
as (this schematic form is modified for negative (3 since in this case fpiu) diverges for 
large u, see (JH])). As Eqs. (|T2l) -(fT^ indicate, the two scaling functions are, to leading 
order in t, identical for any Xi{t) in the range ffTTj) . explaining why the crossover point 
between the two regimes can be chosen anywhere in this range. 

According to the calculation of Appendix A the values of (3, 6 and C depend of 
the initial conditions. We consider initial conditions Go{x) = G{x, 0) which for large 
have an asymptotic form 

Go(|x| > 1) ~ (15) 

Here, a > 1 — b must hold for P{x, 0) to be normalizable. Note that a may be negative. 
If Go{x) decays faster than a power law, we formally take a = oo. A few examples of 
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P(x,0) 


Gq{x) 




a 


A 


5{x - Xo) 


5{x-xq) . _ 
P*(x) 


-1 





-1 


(J^-\x\/xo 


: — : — 

P*(x) 


-1 





-1 




C|x|-(''+i)+^(x) . _ 
P*{x) 


-1 





-1 




c|xr(f-i)+^(x) . _ 


CZ\x\+'^ 


-1 


cz 


P\x) + C\x\-^^+^^ +l{x) 


C|a:|-(''+i)+^(x) 


CZ\x\-^ 


1 


cz 


P*{\x\+xq) + 1{x) 


P*(l^'l+^o)+n^) 1 _ 
P*(x) 


—hxQ\x\~^ 


1 


-bxo 


P*(x)[l + e-l^l/^«] +l{x) 


-\x\/xo , _ 
^ ^ P-W 


g-|a;|/xo 


oo 




(:7P*(x)[l + e-l^l/^°] 


C - 1 + ~ 


1 





C-1 



Table 1. A variety of initial conditions P{x, 0) and the corresponding values of a 
and A according to Eq. ([T5|) . Go (a;) = G(x, 0) is defined by Eq. ([9|), and its leading 
asymptotic behavior for |a;| 3> 1 is presented. £{x) and C are a compactly-supported 
function and a constant whose values change from line to line. They are added to 
ensure the normalization J P{x, 0)dx = 1. The equilibrium distribution P*{x) is given 
in Eq. (0). 



different initial conditions and the corresponding values of a and A are given in Table 
m In Sec. IHwe briefly discuss cases in which Go(x) is not asymptotically symmetric, 
i.e., when the tails of Gq{x) at x — )■ ±oo decay at different rates. 

For this large class of initial conditions, the scaling exponents are given by 

"-"WH f (16) 



and 



S = S(a) = /3{a) + - — - (17) 



For a < 2, the constant C is 

C = —±—^.A, (18) 

while for a > 2, C depends on the full forms of the initial condition and the potential. 
For a = 2 there are logarithmic corrections to Eq. flTUl) . which are presented in Eqs. 
f lXse]) . (1X371) and (IXiOj) of [Appendix A 



We make two comments about the solution (JTOl) . First, we would like to emphasize 
the non-trivial fashion in which the solution depends on the initial condition. The scaling 
functions //? and and the scaling exponent 6 are determined by the value of f3. The 
f3 exponent exhibits a "phase transition" at a = 2, between a regime (a < 2) in which 
(3 depends on the value of a and a regime (a > 2) in which it does not. As discussed in 
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Sec. 12.41 below, it is this threshold phenomenon which ties our scaling solution with the 
problem of velocity selection of propagating fronts. 

The second comment is that this scaling solution is universal, in two ways: it is 
independent of the small-x details of the potential, i.e., of U{x) of Eq. (jlj. It is also 
independent of the small- a; details of the initial condition. When we say below that 
a particular result is universal, we use the term in both these meanings. To be more 
precise, the universal function is ^''^'p^^-^ ^'^^ P{x,t) itself depends on U{x) for small 
values of x, but only through the simple Boltzmann distribution (|H]). It is interesting to 
note that when a > 2, where the solution does not depend on the initial condition, the 
constant C is non-universal, while in the case of a < 2, where the initial condition does 
affect the scaling form, C is universal. 

In the remainder of this section, we motivate these results in order to gain an 
understanding of their origin. To do so, we derive these results in a heuristic fashion 



which, although not rigorous, is more transparent than the calculation of Appendix A 



In this heuristic derivation, we do not presume the results of [Appendix A[ but for a 
single fact: in order to establish the selection mechanism which leads to Eq. f lT6|) (in 
Sec. 12. 4p . we rely on the fact that localized initial conditions Gq{x) (which correspond 
to a = oo) evolve into scaling solutions of the form (fTOj) with /3 = 1. In other words, 
using the scaling solution for localized initial conditions, we are able to find the scaling 
solution for all initial conditions. 

2.3. Scaling solution for \x\ y/t ("large x") 

As we are seeking a scaling form for G{x,t) (Eq. (^) rather than for P{x,t), we start 
by writing down the equation governing the evolution of G. Substituting ([9]) in the 
Fokker-Planck equation ([6]) we find 

dG _ yf(.dG d^G _ _b -^^^ (iQ) 

dt dx dx'^ X dx dx"^ 

Here we have used Eq. ([8]) to deduce that ^^-^^ = —V'{x)P*{x). Our heuristic 
derivation of the scaling solution ffTOj) proceeds by dropping the h{x) term in this 
equation, which is negligible for large values of x. This leads to 

dG{x,t) _ b dG{x,t) d^G{x,t) 
dt X dx dx"^ ' 

which is equivalent to Eq. ([1]). Dropping h{x) is justified below, in Sec. 12. 6[ where we 
establish the universality of the results which we now derive, i.e., their independence of 
the form of h{x). 

The goal of the present subsection is to show that Eq. ( l20l) admits a family of 
scaling solutions. We start by looking for scaling solutions of the form 

G(x,t)=t-^/(^) (21) 
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where the scahng exponent /3 and the function f{u) are to be determined. This 
corresponds to the ansatz 

P(x,t) = P*(x)+P*(x)r'^/(M) + (22) 

for the probabihty distribution. Substituting fl2T]) in the Fokker-Planck equation fl2U]) 
yields a family of ordinary differential equations for /(m), 

r+g-^)/' + /3/ = (23) 
with (3 a free parameter. For every value of /3 this equation has a solution 



(24) 



where iFi is the confluent hypergeometric function 2^, and Ci and C2 are integration 
constants. The three unknown constants /3, Ci and C2 should in principle be determined 
by the two boundary conditions at m = and u = 00 and by the initial condition. 

The study of the small-x scaling solution in Sec. 12.51 below shows that the proper 
boundary condition to consider at m = is 

f{u < 1) ~ v'+K (25) 

Using the asymptotics of the hypergeometric function 2J 

{1 + 0{u^) for M < 1 

^™±^-u-^^ for«»l,r-. ^0,1,2,... , (26) 

^M±^ . «2(-s)e-^ for « » 1, r - s = 0, 1, 2, . . . 

we see that the boundary condition fl25l) implies that C2 = 0, and therefore f{u) = 
C fjiiu), where f/^ is given in ( fT3l) . and we have defined C = Ci. 

Without another condition which may set the values of the two remaining constants, 
C and, more importantly, /9, we are still left with a family of scaling solutions. Note 
that the conservation of probability cannot be used to determine these constants, since 
the scaling ansatz f l2T]) does not hold for small enough values of x. Similarly, the known 
stationary distribution ^ does not provide a boundary condition as all solutions relax 
to it (as can be seen using ( 126|) ). We therefore arrive at the uncommon (although not 
unique, see j25[) situation in which the scaling exponent /3 and the scaling function / 
are determined by the initial condition. This situation confronts us with a problem of 
selection: which of the family of scaling solutions is selected by the initial condition of the 
physical system under consideration? We turn to this question in the next subsection. 



2.4- Stability and the selection of the scaling solution 

In this section we elucidate the selection mechanism which leads to Eq. ( |T6l) . Since 
probability is locally conserved by the diffusion equation (j3]), it is reasonable to expect 
that the relaxation towards equilibrium propagates as a diffusive "front" from the origin 
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towards the tails. If this is so, then at any given time, the tails of G{x^ t) do not yet 
"feel" this front, and they should therefore be given by the initial distribution. By 
matching the tails of Eq. ( 12^ with initial conditions of the form ( |T5l) . the asymptotics 
(!26|) of the hypergeometric function suggest that /3(a) = a/2 and C is given by (ITSll 
when a 7^ 2, 4, 6, . . .. 

According to the exact calculation of Appendix A , the naive argument of the 
previous paragraph is correct only for a < 2. To understand why the argument fails 
when a > 2, we turn to a stability analysis of the scaling solutions and show that those 
with /3 > 1 are unstable to localized perturbations. To this end, we make use of the 
following result which is derived in Appendix A localized initial conditions G{x, 0), such 
as compactly supported ones, evolve at long times to 

G(x,t)~rVi(^) (27) 

where fi{u) is given in f[T^ . This result can heuristically be understood as follows: if the 
initial condition is compactly supported, then it is plausible that the selected solution 
will be the one whose decay at the tails is steepest. The asymptotics (l26ll of the scaling 
solutions show that this is the /3 = 1 scaling function. We note that using the identity 

iF^{A-A-z) = e\ (28) 

one can simplify the expression for the scaling function to fi{u) = -u^^^e""^/^. 

Let us consider a distribution which at some time t is close to a scaling solution of 
the form (l22l) . The distribution cannot be exactly equal to this scaling solution in any 
physical situation: at small enough x's, where the potential deviates from a logarithm, 
the scaling form breaks down. At best, the exact solution is equal to the scaling solution 
plus a small localized disturbance 6P, i.e., 

Pix,t)-P*ix) = CP*{x)t-^fp(^^)+5P{x,t) (29) 

Examining Eq. (l29l) we see that the late-time behavior of the solution will be close 
to the scaling solution fp if the disturbance 6P is negligible compared to it. In other 
words, at late times we can only see scaling solutions which are stable with respect to 
local perturbations. In order to ascertain the stability of the different scaling solutions, 
we should examine how localized perturbations around them evolve in time. Since 
the Fokker-Planck equation is linear, the evolution of such localized perturbations is 
independent of that of the scaling solution. This simplifies the stability analysis: we 
need only to solve the Fokker-Planck equation for localized initial conditions. 

We now use the result (127|) from the exact calculation, and find that at late times, 
Eq. ([29D evolves to 

p*(x)r'^^c(^)+5P(x,t)^p*(x)[cr''/;3(^)+C't-Vi(^)].(30) 

When P < 1 the second term on the rhs is negligible compared to the first, and the 
scaling solution is stable to localized perturbations. On the other hand, when /3 > 1 
the second term dominates the late time behavior. Such scaling solutions are unstable 
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and can never be observed in physical systems. We thus see that initial conditions of 
the form (|T5|) with a < 2, which "excite" scaling solutions with /3 < 1, evolve to a 
scaling solution which depends on a. "Most" initial conditions, however, evolve towards 
the marginally stable scaling solution of /3 = 1. By "most" we mean that the basin 
of attraction of the /3 = 1 solution in the space of all initial conditions has a higher 
dimension than the basins of attraction of solutions with any /3 < 1. Note that when 
a > 2, the constant C is determined by the localized perturbation rather than the tail 
of the initial condition, and therefore it is not given by Eq. (in 



2.5. Scaling solution for \x\ <^ \lt ("small x") 

In this section we show that when |x| ^ the probability P{x,t) is also given in the 
long time limit by a scaling form. This scaling form is different from the one discussed 
above, but it, too, depends on the initial condition. Surprisingly, this scaling form is 
universal: it depends on the full details of the potential V{x) only through the stationary 
distribution (|8]) which multiplies the scaling function. 

Unlike the large-x scaling form, the small-a; scaling form depends on the boundary 
condition at the origin. As mentioned above, we assume in this section that the 
probability current at the origin (defined in Eq. (j3])) vanishes at all times: J{0,t) = 0. 
This boundary condition translates into 
dG 



dx 



= (31) 

x=0,t 

as long as P*(0) 7^ 0. The results for other boundary conditions are discussed in Sec. HI 
From the calculation of Sec. 12.31 we already know that at a;'s which scale as y/t, 
G{x,t) is given by fl2T]) and f lT3|) with /3 which depends on the initial condition. We now 
examine the solution at x which scale as t''^, with < 7 < |. To this end, we look for 
scaling solutions of the form 



X 



G{x,t)=t-''^g,i^-), (32) 

which we call "the solution at scale f^" . 

We begin by considering the unsealed solution G{x, t) itself (this is the case 7 = 0). 
Substituting the ansatz 

Gix,t)=t-'°goix), (33) 
in the Fokker-Planck equation f|T9|) yields 

g'^ix) - V'{x)g'o{x) = -6ot-'go{x). (34) 

For t ^ 1 the term on the rhs becomes negligible^ We thus arrive at the simple equation 

I More precisely, we expand G{x,t) in a power series in t^^: G{x,t) — G3t^^[aQ{x) + ai{x)t^^ + 
a2{x)t^'^ + ■■•], which we substitute in ([T9l) and solve separately at each order. As is shown below, at 
the zeroth order we find ao(a;) = 1. The next order gives ai{x) = C' — S dy dzexp{V{y) — V{z)), 
which for a; ^ 1 is approximately ai{x) ~ x'''^^. This means that as long as a; <C t^/^^+i) the 
approximation go{x) — C3 is valid. 
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(7q(x) — V(x)5'q(x) = which can be integrated, yielding 

g,{x) = C^ + C^ r e^'^y^dy, (35) 



where C3 and C4 are integration constants. The boundary condition (13T!) imphes that 
fl'o(O) ~ 0' which means that C4 = 0. Therefore, for values of x which are small enough, 
G{x,t) = C3t-^°. 

We now proceed to examine the solution at scales with < 7 < |. At late times, 
X ^ ^ 1, and we can replace V{x) ~ blogx. Substituting the ansatz fl52]) in the 
equation (1201) yields the ordinary differential equation in the scaling variable z = xt~'^ 

g'!^{z) - -J^{z) = -bzg'^iz) + S,g,iz)r^'-'^\ (36) 

As before, we assume that the two terms which are proportional to t^^^^^i) q^-q negligible 
at large times and we drop them. The validity of this assumption will be examined below. 
We are left with the equation: 

g'^iz) - -J^iz) = 0, (37) 
whose solution is given by 

g,{z) = C, + C,z''-\ (38) 

The picture that emerges is that at every scale the solution is either a constant C5 or 
a power law Cqz'^'^^. At a single intermediate scale, both C5 7^ and Cq 7^ 0. Continuity 
at small x implies that if C5 7^ then C5 = C3. With an abuse of notation, we shall from 
now on denote the exponent of this special intermediate scale by 7. At this intermediate 
scale, we expect the solution for large values of z to coincide with the small u behavior 
of the solution in the x ~ a/^ region (see Eqs. f|T0|) - f|T6|) ). This yields the condition f l25|) . 
Using G{ut^/^,t) ^ Ct-(^^^\^+^ (see ([HD), we find that Cg = C, and 

5^(a) = /3(a) + (6+l)(i-7). (39) 

We note that for the scaling solution fl38|) . the two terms that were neglected when 
passing from Eq. fl5^ to (1571) are indeed negligible as long as 2; -C t^'"/^ or equivalently 

X < 

We are now left with the problem of ascertaining the values of the two remaining 
undetermined constants 7 and C3. These can be found with the help of the conservation 
of probability (we once again make use of the boundary condition (!3TI) ). Choosing 

t"' ^ Xi(t) ^ t^/^, we can write 

pco poo 

0= [P{x,t) - P*{x)]dx = P*{x)G{x,t)dx = h{t) + hit), (40) 
Jo Jo 

where we have defined 

/•xi{t) rxi(t) 

h{t)^ P*{x)G{x,t)dx^t-^-' P*{x)[Cs + Gx''+H~^^''+^^]dx, 

Jo Jo 



00 



h{t)= j P*{x)G{x,t)dx^t~^ I ^x-^fJ^)dx. (41) 
Ai(i) Jxi{t) Z Vvt/ 
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In these equations we have substituted the small-x and large-x scahng solutions for G. 
To leading order in t, Eq. (gO]) gives C^t'^-^ = -f§ u'^ fp{u)dut-^-^^-^'^/'^ , where 

2 P*i )d < ^ for diffusion on positive half-line 
Jq I 1 for diffusion on entire real line 



Using ( 139|1 we deduce that 

T 



7 = (43) 



and 



2C p , ,^ C 4(6 + 1) 

= io ^ ^^(^^'^ = ■ 2/3(a) + 6 -l - 



To sum up, we see that for any Xi{t) in the range (fTTjl . 

G{\x\ < xi(t),t) ^ J^), (45) 

where 

9^i^) = -Tz§lh) + ^''^'- (46) 
Once again we find a scaling solution which depends on (the tails of) the initial condition. 
On the other hand, just like the scaling solution at large-x, this solution is essentially 
independent of the full details of the potential V{x), which only serves to determine the 
stationary solution P*{x) and, when a > 2, the constant C. 

We emphasize that the analysis presented above holds for all |x| < Xi(t), including 
the region around the origin where the potential is not logarithmic. For any fixed x 
(which does not scale with t), Eqs. fH5|) and fH6|) agree with rigorous results obtained 



for discrete random walks in a logarithmic potential 26|, |27 



2.6. Universality of late-time scaling solutions 

In this section, we establish the universality of the results of Sec. 12.31 That is, we show 
that they depend only on the logarithmic tail of the potential, and not on the h{x) 
correction term of Eqs. (JTj) and (|T9l) . Moreover, our argument demonstrates that the 
details of the initial condition near the origin are also irrelevant for the large-x scaling 
form. 

To establish the required universality, we rescale x, t and G(x, t) by defining a 
rescaled function 

Gp^{x,t)=K^^G{Kx,KH). (47) 

Thus, up to a normalization factor which depends on /3, the rescaled function G\{x,t) 
is equal to G at time K^t as seen at the spatial scale Ax. The equation for the evolution 
of G\ may be straightforwardly obtained by substituting the definition ( H71) into Eq. 
(dS]), yielding 

dGK{x,t) b dGA{x,t) d^G^{x,t) 

-^r- = -x^^ + ^^(^^^^^ + dx- ' ^^^^ 
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where h\{x) = h{Ax). 

The solution G{x, t) at a given late time t S> 1 can be obtained in two ways: either 
by propagating the initial condition according to Eq. (fT9|) . or by rescaling the initial 
condition, propagating it according to Eq. (148|) to (rescaled) time 1 and rescaling back. 
In the second way, the correction h{x) is negligibly small. According to this procedure, 
we obtain the scaling limit by replacing A with \/t: 

t^G{ut^'\ t) = 1) ^ Goo(n, 1) = f{u). (49) 

The exponent /3 must be chosen appropriately so that the t — )■ oo limit exists and is not 
zero. The limiting function Goo{x,t) evolves according to fj48|) with 

h^{x) = lim h^ix) ~ lim (Ax)"'" = (50) 

A— >-oo A— >-oo 

for any x 7^ 0, see ([Tj). Therefore, the scaling limit of the original Fokker-Planck equation 
does not depend on h{x). 

The initial condition for the rescaled problem fHS]) is 

G^{x,0) = lim Ga{x,0) = lim A'^^G{Ax,0). (51) 

A— s>oo A— s>oo 

This limiting initial condition Goo{x,0) is in many cases a singular function, similar to 
the Dirac (5-function but with a different type of singularity. If the tails of G{x, 0) at 
X — )• ±00 decay with \x\ faster than algebraically (e.g., exponentially), then Goo{x,0) is 
zero when x ^ and is singular at the origin|§| The exact details of the initial condition 
around the origin are lost in the limit which yields Goo{x,0). The tails of the initial 
condition may, however, affect Goo: an initial condition which decays algebraically as in 
Eq. (fT5|) is rescaled to Ga(x, 0) ~ A'^'^~"-A\x\~"'. If /3 = a/2, the rescaled initial condition 
Goo has the same algebraic decay as G. We see that initial conditions may affect the 
scaling solution only through their tails, and that for initial conditions with power-law 
tails, a limit of flSTl) exists only when choosing P < a/2 (compare with (|T6|) ). 

The rescaling argument which we have presented in this section is inspired by the 



renormalization group (RG) techniques used by Goldenfeld et al. [28| and by Bricmont 



and Kupiainen |29| to analyze nonlinear partial differential equations. Here we have used 
their method to analyze a linear equation which is inhomogeneous in space. From the RG 
perspective, the rescaling transformation fH7|) can be viewed as an RG transformation 
that has a one-parameter family of fixed points /^(m). The scaling limit of the original 
equation is determined by the appropriate fixed point, which is not affected by the 
addition of h{x). Therefore, the h{x) term in the equation is irrelevant and the scaling 
solution is universal in the RG sense. 

§ Consider for example a symmetric localized initial condition G{x, 0) which is negative at the origin, 
becomes positive at = 1, and is exactly zero for \x\ > 2. In this case G\{x,0) is somewhat similar 
to S"{x), the second derivative of the Dirac S function, but might be either more or less singular than 
S"{x): if (j){x) is a smooth test function, then / (j){x)GA{x,0) exp[~VA{x)]dx ^ 0"(O)A2^+''~^, which in 
the limit A — > 00 might diverge or vanish, depending on the sign of 2/3 + 6 — 3. Here, Va is the rescaled 
potential, defined by V^(a;) = |(1 + hA{x)). 
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3. Comments on the scaling solution 

In this section we comment on the scaling solution derived above and discuss it in some 
broader contexts. 

3.1. The scaling solution and incomplete self-similarity 

The dependence of the scaling exponent /3 on the initial condition signals the failure of 
dimensional analysis. The latter is easily seen to predict incorrectly that (3 = The 
reason for this failure is the following. The prediction of dimensional analysis for the 
diffusion equation rests crucially on the conservation of probability [25]. When 6 > 1, 
the limiting rescaled equation ([T]) does not conserve probability at the origin because 



of the singularity of the potential there 2l|, |22(|. It is the U{x) term in the potential 



dl]) which guarantees conservation of probability, and rescaling it away (as was done in 
Sec. 12. 6p yields a singular limit. In practice, this means that at any finite time t, no 
matter how late, corrections due to the potential U{x) inevitably affect the form of the 
solution at small enough values of x. The scaling form (1221) . in which these corrections 
are not taken into account, should not be expected to conserve probability by itself, 
and therefore, the scaling exponent /3 cannot be found by dimensional analysis. In the 
terminology of Barenblatt, the scaling solution to our problem exhibits self-similarity of 



the second kind (see 25l |) 



It is interesting to note that when b < 1, the solution of the diffusion equation ([2])- 
(|3]) approaches a self-similar solution of the first kind, i.e., one whose scaling exponents 
can be determined by dimensional analysis. The scaling solution in this case was found 



m 



19l . |20| . Since there is no equilibrium distribution when 6 < 1, one cannot define 



G{x, t) according to Eq. iQ. Nonetheless, one may look for scaling solutions of the form 
P{x,t)r^x-h-^f[^). (52) 

Dimensional analysis (or, equivalently, the conservation of probability) dictates as before 
that /3 = which in this case is indeed the correct value, regardless of the initial 
condition. For example, when 6 = 0, i.e., in the simple case of free diffusion, one obtains 
the well known result [i = \. 

3.2. Selection and propagating fronts 

A selection mechanism similar to the one described in Sec. 12. 4[ by which most initial 
conditions evolve into a marginally-stable state, is well known to exist in several other 



problems [17|, l30N32| . Many of these problems can be expressed as propagation of fronts 
into unstable states [l3]. A well-studied example is given by the non-linear diffusion 



equation which was studied originally by Kolmogorov, Petrovsky and Piskunov 33| and 
by Fisher 34| : 



dt dx^ 



+ (53) 
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Their original works concern the spreading in space of an advantageous mutation in a 
population. In this context < t) < 1 describes the fraction of individuals located 
at point X who posses an advantageous gene. This equation admits two stationary 
homogeneous solutions: an unstable solution = and a stable solution 0(x) = 1. 
Any localized initial perturbation around the = state grows into two traveling waves 
propagating outwards with an asymptotically constant velocity. This velocity of front 
propagation cannot, however, be easily determined, as Eq. (|53|) has a traveling wave 
solutions 4>{x,t) = fv{x — vt) for every possible velocity v. 

The selection mechanism for the problem of propagating fronts has strong 
similarities to our problem of diffusion in a logarithmic potential. It is possible to 
show [l3] that, similarly to our problem, the selected front solution depends on the tails 
of the initial condition: if 0) ~ e~'^^, then the asymptotic velocity is f (A) = A + 1/A 
for A < A* = 1, and is v{\*) = 2 independent of A for steep enough initial conditions, 
i.e., when A > A* (the latter case includes localized initial conditions, i.e., those which 
have a compact support). Moreover, all traveling wave solutions with v < v{X*) are 
unstable to small, localized disturbances. Notice also that both in our problem and 
in the problem of front propagation, stable solutions decay monotonically at the tails, 
while unstable solutions decay at the tails through oscillations (in our problera, this is a 



property of the hypergeometric function (fT3|) : for propagating fronts see, e.g., [17|]). The 
marginally stable solution, into which localized initial conditions evolve, is the solution 
with the steepest tail which is still monotonic. 

The similarity between our problem and the selection of propagating fronts is 
furthered by noticing that, by a simple change of variables, scaling solutions in general 
can be thought as traveling waves |25|]: by defining ^ = logo; and r = logt, any scaling 
solution can be expressed aijj] 

^~^f{^) = e-^^fie^~n = e"^-0(e - vr), (54) 

which is a traveling wave solution (whose overall height might shrink or expand with 
time, depending on the sign of f3). Note also that a power law tail of the initial conditions 
f lT5|) implies an exponential tail in the traveling wave variables: Ax~°' = Ae~°'^. 

A few peculiarities of the selection problem posed by Eq. (E]) should be mentioned. 
First, unlike in the problem of propagating fronts, the velocity of the traveling wave 
( 154|) which corresponds to the fast front solution ( ITO!) and (fT6l) is independent of the 
selected solution: it is always v = 1/2. Instead, it is the exponent (3 which is selected 
by the initial condition. In addition, as discussed above, localized initial distributions 
for the diffusion equation correspond to initial conditions ( |T5l) with a = 0, which, in the 
context of selection, are not localized (in other words, when the distribution P{x, 0) is 
localized, then G{x, 0) is not localized, and /3 = is selected rather than the marginal 
value 1). This is in contrast with the many problems of selection in which the generic 
initial conditions which are natural to consider are the localized ones. Furthermore, 

II This transformation is only valid for x > 0. One can separately transform the negative x scaling 
form into a traveling wave solution by defining ^' = \og{—x). 
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other "non-steep" initial conditions are physically relevant in many situations, as will 
be discussed in Sec. [51 In other words, unlike many other selection problems, Eq. ([3]) 
naturally leads us to study those cases in which the solution does depend on the initial 
condition (another such exception is found in [35]). 

Another difference of the diffusion problem from most known problems of selection 
lies in the fact that the scaling solution of the diffusion problem is made up of two scaling 
functions. As mentioned above, the scaling form of the slower front is determined by 
that of the faster one, and hence it is also selected by the initial condition. In the 
language of traveling waves, this corresponds to a case in which two moving fronts exist, 
propagating at different velocities. While there are systems which are known to develop 
two fronts selected by a marginal stability mechanism, we are not aware of a case in 
which the velocities of the two fronts are related to each other by an expression akin to 
Eq. (HZ]). 

An interesting feature of Eq. ([3]) is that, unlike other problems where selection takes 
place, this equation is linear, yet not homogeneous in space. The linearity of Eq. ([3]) 



enables the derivation of an exact solution (as is done in Appendix A ), and thus assists 
in analyzing the selection mechanism in detail. It should be noted that while many 
problems of selection are conjectured to be governed by a marginal stability criterion, 
a rigorous proof of this fact is rarely known. The simpler linear example provided by 
Eq. ([3]) might help to shed light on the common mathematical structure governing these 
similar problems. 

3. 3. Relation with previous results 



We briefly comment on the relation of our results to those of |18|-|20| . There, a scaling 
solution of the form (152|) rather than (12T|) was sought. For such a scaling solution, /3 
must be equal to zero and 

/(m < 1) = 1/Z + 0(m) (55) 

must hold, since P{x,t) should eventually converge to the steady state distribution ([H]). 
The scaling function / satisfies the same differential equation ( l23i) as /, whose solution is 
(^^. Using the asymptotics of the hypergeometric function (l26ll . the boundary condition 



(1551) together with f{u — )■ oo) — determine the constants Ci and C2, and the scaling 
solution is found to be f{u) = F(^, ^)/ZF(^), where T(a,x) is the incomplete 
F-function. 

Examining the general solution (fTOj) . ( [T3|) . (fT6|) . and (fT8|) . and using properties of 



the hypergeometric functions 2^, it can be verified that when a = and A = —1, the 



scaling exponent (3 indeed equals zero and the large-x scaling solution reduces to the 



results of |18H20|. This case includes the large class of initial conditions P(x, 0) which 
decay to zero faster than a power law, e.g., P(x,0) = S(x — Xq) (see Tabled]). Other 
initial conditions, however, select different values of (3 and lead to a scaling function 
different from the one considered previously. 
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4. Non-conserving boundary conditions 

So far, we have concentrated on solutions of Eq. ([3]) with no-fiux boundary conditions 
at the origin, i.e., we have assumed that the current of probabihty at the origin J(0,t) 
is zero at all times. In this section, we describe what happens for J(0, t) 7^ 0. Such a 
situation arises in two different scenarios: (i) the distribution P{x,t) is defined only 
for X > and the boundary condition at the origin allows J(0, t) 7^ 0; (ii) x is 
unbounded, but the initial condition is asymptotically non- symmetric, i.e., G{x — )■ 
±00,0) ~ A^x'""^ with a+ 7^ a_ or 7^ A^. We focus here on the first scenario, 
and only briefly describe what happens in the second. 

For concreteness, we discuss a specific choice of boundary condition at the origin: 
an absorbing boundary condition, i.e. P(0,t) = 0. This boundary condition arises 
naturally in many physical problems, especially when studying first-passage properties 
of the dynamics (see Sec. 15. 3p . Other boundary conditions (e.g., P{0,t) = Pq where Pq 
is a constant) can be treated in a similar manner. We remark that diffusion on the half 
line X > with an absorbing boundary at the origin is equivalent to diffusion on the 
entire real line with an initial condition which is antisymmetric. This suggests that the 
case of an absorbing boundary can be treated similarly to the unbounded x which we 
have considered in previous sections. We do not follow this alternative route below, as 
we seek a derivation which can easily be generalized to other boundary conditions. 

When probability is not conserved at the origin, the eventual steady state which the 
system reaches need not be P*{x) (which by definition ([8]) is normalized to 1). Thus, 
we need to redefine G{x,t), as we expect to find a scaling form for solutions which 
eventually relax to zero. We therefore define 

P^{x) = lim P{x,t) (56) 

t— >oo 

to be the steady state which the system eventually reaches, and generalize the definition 
of G to 

G(x,t) = ^^;I^ (57) 

(compare with (^). Note that Poo{x) depends on the boundary condition at the origin. 
For instance, for a reflecting boundary condition Poo{x) = P*{x), while an absorbing 
boundary results in Poo{x) = 0. The definition (|57l) allows us to consider both cases on 
the same footing. 

The parameter a is defined by the tails of G{x,0), which depends by definition 
on the boundary condition (see Eq. dSZ])). Therefore, the same initial condition 
P{x, 0) may result in two different values of a when considering two different boundary 
conditions (conversely, one may say that for different boundary conditions, the same 
initial condition G{x,0) corresponds to different initial distributions P{x,0)). A few 
examples of different initial conditions and the corresponding values of a and A in the 
case of an absorbing boundary are presented in Table O 

Examining the argument of Sections 12. 31 and [2l4t we see that the boundary condition 
at the origin does not play any role in the derivation of scaling form at large values 
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P(x,0) 


G{x 


,0) 




a 


A 


5{X - Xq) 


5{x-xo) 
P*{x) 







oo 




(J^-\x\/xo 


P*{x) 


r-u 




oo 






C|a;|-(''+i)+£(z) 


rsj 




1 


CZ 




C|a:|-(''-i)+^(x) 
P*(:r) 


rsj 


CZ\x\+^ 


-1 


CZ 




. 1 C|x|-(''+i)+^(x) 


rsj 


1 





1 


P*{\x\+Xq)+1{x) 


P*{\x\+XQ)+t{x) 
P-ix) 


rsj 


1 





1 


P*(x)[l + e-l^l/^'o] +^{x) 


1 _L p-|a;|/a;o 1 


rsj 


1 





1 


CP*(x)[l + e-l^l/^«] 


C + Ce-l^l/^° 




C 





c 



Table 2. A variety of initial conditions P{x, 0) and the corresponding values of a 
and A according to Eq. ()15p for a system with an absorbing boundary at the origin 
(in which case Poo{x) — 0). Go{x) is defined by Eq. (|57|) . and its leading asymptotic 
behavior for |a;| ^ 1 is presented. £{x) and C are a compactly-supported function 
and a constant whose values change from line to line. They are added to ensure the 
normalization J P{x, 0)dx = 1. The equilibrium distribution P*{x) is given in Eq. ([5]). 
In some cases, the values of a and A for the same initial condition might differ when 
the boundary condition is changed (compare with Tabled]). 



of X. We can therefore conclude that the large-x scahng form is independent of the 
boundary condition at the origin. This conclusion is supported by the exact calculation 



of Appendix A 



The small-x scaling function, on the other hand, does depend on the boundary 



condition. As in Sec. 12.51 we start by considering an ansatz (1331) for the unsealed solution 
G{x,t), and obtain Eq. f p5|) . When the origin is absorbing, the boundary condition on 
G is G{x, 0) = 0, from which we deduce that C3 = 0. We therefore have 

G{x<xi{t),t)^Ct~'gi,{x), (58) 

where 

~gp{x) = C, r e^^y^dy, (59) 



^0 

and Xi{t) is in the range 1 ^ Xi{t) ^ \/t (compare with fH5|) - fH6l) ). The constant 
G is the same as in Eq. flTS]) . and = G^/G. From f l5I?]) together with the form ([2]) 
of the potential it is seen that ^^(x) ~ x^^^/{b + 1) for x ^ 1. Matching the large-a; 
asymptotics of (158|) with the small-w asymptotics of (fT3l) yields 

6{a) = /3(a) + ^ (60) 

and C*4 = 6 + 1 . Note that the new form of the small-x scaling function (1591) is no longer 
independent of the small-x details of the potential V{x). 
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Figure 2. A schematic representation of the solution P{x,t) (Eq. (jBTj) ) at a given 
late time t 3> 1 (not drawn to scale) for a system with an absorbing boundary at the 
origin. The red hollow line represents gp{x)t~^ , the solution at small values of x, while 
the blue solid line represents the large-x scaling form fi3{x/t^^^)t~^. The boundary 
condition may modify the value of a, but it does not have any further effect on the the 
large-x scaling form (compare with Fig. [T]). 

To sum up, changing the boundary condition at the origin affects the scahng 
solution in two ways. First, it entails a change in the definition of G{x, 0), which might 
alter the value of a. Second, it modifies the small- a; scaling function. Importantly, the 
large-x scaling function and the scaling exponent /3 remain unchanged. For the case of 
an absorbing boundary, these changes are summed up in the final scaling form of the 
solution 

f ^^(a)(x)r^» for X <xi(t) 
P{x,t)^PUx) + CP*ix)- I , (61) 

[ //3(a)(^i72)r^("^ forx>Xi(t) 

where, g is given in Eq. f l59|) and 6{a) by fl60|) . This solution is depicted schematically 
in Fig. [2 

Finally, we briefly comment on diffusion on the entire real axis with non-symmetric 
initial conditions G{x — )■ ±oo, 0) ~ A±x~"'^ (scenario (ii) above). In this case, there are 
two different "fast fronts" which propagate from the origin to ±oo: a scaling function 
for X ~ y/i and another for x ~ —y/t. Each of these is selected by the corresponding 
tail of the initial condition. Similarly, there are two "slow fronts" (i.e., small-x scaling 
functions), each one overlapping with the corresponding large-x scaling function. A 
calculation similar to that of Sec. 12.51 can be repeated, leading to solutions of the form 
fl32|) . fl38l) and f p9|) . Four unknown variables remain: 7-1- and C5-1- (the values of 7 
and C5 for the positive-x and negative-x scaling functions). These can in principle be 
determined from two equations: the conservation of total probability, and continuity of 
the probability current at a; = 0. 
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5. Applications 

In this section, we present applications of the new theoretical results which have been 
derived above. In particular, we give examples of several problems in which the 
dependence of the scaling form on the initial conditions plays an important role. 

As discussed above, when probability is conserved, a large class of initial 
distributions (including localized ones) correspond to a value of a = (see Eqs. ([9]) 
and ( fTSll and Table [1]). For these initial conditions the distribution evolves to the /3 = 



scaling form, which is the one previously obtained in [18H20[. An inspection of Eqs. 
([9]) and ( |T5|) and of Table [H reveals that initial conditions with a 7^ can be divided 
into two broad classes: for a < 0, the exponent a yields the leading decay of the tail 
of the initial distribution. On the other hand, when a is positive, the tail of the initial 
distribution approaches the equilibrium distribution P*(x); in this case, the leading 
decay of the initial distribution as x — ?■ ±00 is that of P*, and a determines the sub- 
leading correction to P* . Below we consider examples of both classes. In Sec. 15. we 
describe an experimental protocol by which initial conditions with negative values of a 
can be obtained, and we propose a cold-atoms experiment which, using this protocol, 
could measure the predicted dependence of the relaxation on the initial condition. Initial 
conditions belonging to the second class may at first sight seem unnatural in physical 
circumstances, as they require fine-tuning the initial distribution. In Sec. 15.21 we show 
that this is not necessarily the case, and explain how initial conditions with a = 1 arise 
naturally in the calculation of current correlations in the zero-range process, a stochastic 
model of particle transport. 

When the boundary condition at the origin is absorbing, on the other hand, the 
value of a is always determined by the leading decay of the tail, no matter what the 
initial distribution is (see Eq. ( [57|) and Table [2]). Therefore, no value of a requires fine 
tuning of the initial distribution. In Sec. 15.31 we provide one example of such a system: 
we explain why the dynamics of loops in a denaturating DNA molecule is described by 
Eq. ([T]) with an absorbing boundary, and show the implications of the dependence on 
initial conditions to the analysis of results of single-molecule experiments. 

For the sake of completeness, we provide in Sec. 15.41 a review of some other systems 
which are described by Eq. ([T]), to which our results may be relevant. 



5.1. Initial conditions with a < and atoms in optical lattices 

Equations ([9]) and ( fTSl) indicate that the tail of the initial distribution for a < is of 
the form 

P(x, 0) ~ Ax'^" with < 6. (62) 

Here, a = fi — b. Thus, such initial distributions can be relatively easily generated in 
physical situations. This observation straightforwardly suggests a protocol by which one 
can observe the dependence of the relaxation dynamics on initial conditions. For the 
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sake of concreteness, we present this protocol in the context of cold atoms trapped in 
optical lattices, where the dependence on initial conditions can be tested experimentally. 

When cold atoms are placed in optical lattices, their momentum performs a diffusion 
which, in the semi-classical regime, is of the form ([1]) where x represents the momentum 
[l|. In recent years, this momentum diffusion has received both theoretical and 
experimental attention due to the power-law distribution and "anomalous" dynamics 
to which it gives rise j^^, 36 1. The Fokker-Planck equation for the semi-classical 



probability distribution W{p,t) of an atom with momentum p at time t is 

dW{p,t) d r , dW{p,t) ^ 

Qt = ^ [-FipWip, t) + D{p) J , (63) 

where, in appropriate units, -F(p) = — = ~^ + 0{p^'^) is the cooling "friction" force. 



and D{p) = 1 + is a momentum-dependent diffusion coefficient The parameters 
b and D are determined by the depth of the optical lattice, which may be controlled in 
an experiment by the intensity and detuning of the optical lattice. When D ^ 1 Eq. 
( 163|) is of the form (|2])-([3]). The equation can be brought to this form even when D is not 
negligible, by the standard transformation q{p) = D/D{p')dp' = p + Darctan(p) 



231]. The transformed equation reads 



where once again F{q) = — - + 0{q~^). For convenience of notation, we will assume 
below that D 1 and study Eq. ( 163|1 . 

The experimental protocol to observe the "anomalous" scaling suggested by 
equations f|T0|) - f|T6|) is rather straightforward. For any given value of the parameter b, the 
stationary distribution of momentum is given by W^{p) = Z^^e^^^^P"* = Z^^{1 +p^)~^/^ 
(where we have made the dependence on the parameter b explicit in our notation). 
In an experiment, the parameter b can be controlled by changing the depth of the 
optical potential. The following two-step procedure would generate an appropriate 
initial condition with negative a: (1) a state with momentum distribution W{p, 0) = 
Z^^ai^ + p2)~(^+")/2 with some 6 > 1 and 1 — 6<a<0is prepared by setting the 
parameters of the experiment to a value which corresponds to 6 -|- a, and allowing the 
system to equilibrate; then (2) at time t = the parameters are rapidly changed from 
b + a to b. Following this "quench", the distribution W{p) or one of its moments is 
measured as a function of time. For instance, if 6 + a > 3, one may measure the 
variance of the momentum (p^) (which is proportional to the mean kinetic energy of the 
atom), which is predicted to decay as 

{p'it)) - {p\00)) = I dpp'[Pip,t) - P^ip)] ~ I u'-'fa/2iu)du. 

(65) 

Although we have presented this experimental protocol in the context of cold atom 
experiments, it could be used in many other physical contexts as well. 
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We remark that unlike many cases in which fat-tail distributions lead to an 
anomalous time-evolution, in the case which we discuss here there is no requirement 
that any particular moment of the initial distribution diverges. In fact, any particular 
moment of the initial (or final) distribution can be guaranteed to be finite by selecting 
h large enough with a fixed value of a. 

5.2. Initial conditions with a = 1 and current correlations in a critical zero-range 
process 

Another way to generate initial conditions with a 7^ without fine-tuning the parameters 
of the initial state is to prepare the system initially in a translate of the equilibrium 
distribution, i.e., P(x,0) = P*{x + Ax) for some Ax. In this case, the initial condition 
corresponds to a = 1 (see Tabled]). Such a situation may be realized experimentally if 
it is possible to displace the confining logarithmic potential. 

In this section we present a different case in which such an initial condition 
arises. The problem we shall address here is the calculation of stationary two-time 
correlations of particle currents in a zero-range process (ZRP), a stochastic model of 
particle transport exhibiting real-space condensation. 

In the ZRP which we consider, N = pL particles hop on a one- dimensional lattice 
of L sites with periodic boundary conditions (p is the density of particles). The particles 
can only move in one direction. The defining property of the model is that the rate of a 
jump from site i to i -|- 1 is a function only of the number of particles in the departure 
site. We denote this rate by w{ni). This non-equilibrium model of interacting particles 
has been studied extensively in recent years. For certain choices of the hopping rates, the 
model exhibits a condensation transition whereby, when the density is increased above 
a critical density pc, a finite fraction of all particles resides in a single site (selected 
at random). For reviews of this condensation transition and other applications of the 



model see 137-39 



We consider a ZRP at the critical density, and examine correlations of the current 
flowing across a single site. We concentrate on hopping rates which for large n have the 
form 

uj(n) = l + - + 0(n'^). (66) 
n 



These commonly studied rates give rise to condensation when b > 2 [40|. In the 
thermodynamic limit (when L — )■ 00), the arrival of particles into any site is a Poisson 
process with rate 1 which is independent of the the process of particles departing from 



the sitq^ll [4l|- The occupation probability of the site P{n) evolves according to the 
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master equation 

d 

-^P{n) = P(n-l) + w{n+l)P{n+l) - [1 + w{n)]P{n) ^ 
^ In a system of finite size L 3> i, the arrival process is approximately Poisson on time scales t <^ L 



4l|. Therefore, the result obtained below (Eq. [73]) is correct for finite systems in the intermediate 



asymptotics regime oi 1 <^ t <^ L. 
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.|_[^(,,0,„-))pH.^, (67, 

which is of the form of Eq. (EIl^dZD- It is straightforward to verify that the steady-state 
distribution is 

^•(») = in^ = 5"-'(l + 0(„-)) (68) 

k=l ^ ' 

where Z is & normahzation constant and Z is non-universal and depends on the full form 
of the rates w{n). For w{n) = 1 + 6/n, for example, it can be shown that Z = h/{b — 1) 
and Z = [{b - l)r{b)]-\ 

Having presented the model, we now present the specific problem which we wish to 
study, and show how the results of previous sections can be used to solve it. Our task 
is to calculate the correlation function 

C(t) = a„,out(t) = (ji„(0)jout(t)) - t > 0, (69) 

where jij^{t)dt is the number of particles arriving at the site between time t and t + dt, 
jout(t)dt is the number of particles departing from the site during this time period, 
and j = {jin{t)) = {joutit)) = ^P*(n)u'(n) = 1 is the mean current in the steady 
state. Angular brackets denote an average in the steady state. We may similarly define 
the correlation functions Cin,in(i^), C'out,out(i) and Cout,in(^), but these are all equal to 
zero: both the arrival process of particles entering the site and the departure process of 
particles leaving it are Poisson processes0 and the arrival process is independent of the 
departure process. To simplify notation we shall from now on drop the subscripts and 
denote C{t) = Cir,,ont{t). 

Although the exact steady-state distribution of the model can be calculated for 
any jump rates, little is known about two-time correlation functions such as C(t), even 
in the steady-state. We now show that the long time asymptotics of this correlation 
function can be found using the scaling solution f lTOj) and f|T6|) of Eq. fl67|) with a = 1. 
To this end, we note that (jin(O)jout(^)) is given by a product of the rate with which 
a particle enters the site at time 0, (which is 1) and the conditional rate with which 
a particle leaves the site at time t given that a particle has entered at time zero. The 
latter rate depends on the (conditional) occupation of the site at time t, and therefore 
the correlation function is 



C(t) = ^P*(n) ■ 1 ■ P{m,t\n + 1,0) - P 



[m] 



wim) 



oo 

'^^P{m,t\Po) - P*{m) w{m), (70) 



m=l 



where P{m,t\n,0) is the conditional probability to have m particles in the site at time 
t given that there were n at time 0, and in the last equality we have introduced the 

+ It is not a trivial statement that the departure process is a Poisson process. In the field of queueing 
theory, this statement is known as Burke's theorem, see 42l |. 
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notation 

oo 

P{m,t\PQ)~^PQ{n)P{m,t\n,{)) with Pq^ti) = P* {n - I) . (71) 

For large n, the initial condition Po(^) satisfies 

P^^n) = P*{n)\l + bn-^ + 0{n-^)\, (72) 



and therefore it is of the form ( |T5l) with a = 1 and A = b. 

The calculation of C (t) now proceeds by substituting the appropriate solution (ITOl) 



in Eq. (1701) and evaluating the sum. We carry out this calculation in Appendix B This 
calculation turns out to be somewhat subtle, as the leading terms in t exactly cancel 
out, and the decay of correlations is determined by the next-to-leading term. We note 
here that the cancelation of the leading-order terms can only be established using both 
the small-x and large-x asymptotic regimes. The result of the calculation is 



_ 6+1 



where Z is defined in Eq. f l68|) . 

5.3. Absorbing boundary conditions and dynamics of denatured DNA loops at criticality 

The analysis of Sec. S] has revealed that initial conditions with any value of a can be 
achieved without fine-tuning when the boundary at the origin is absorbing (see Table 
12]) • Absorbing boundary conditions arise naturally when studying first-passage problems 
such as the mean time it takes a diffusing particle to reach the origin from a given initial 
condition (see for example 43|). In this section, we discuss one such example in the 



experimental context of the dynamics of denaturing DNA molecules. 

It is well known that when the double stranded DNA molecule is heated, it 
undergoes a denaturation phase transition in which it separates into two single strands. 
The nature of this phase transition has been debated over the years. Many of the 
theoretical studies of this transition are based on the model of Poland and Scheraga 



44j-|46| (for recent reviews see [471]). These studies model the DNA molecule as an 
alternating sequence of bound segments and denatured loops, or bubbles. The bound 
segments are considered rigid, with each bound pair contributing a negative energy of 
— e in the case of homopolymers, while the shape of open loops may fiuctuate and thus 
contribute to the entropy of the molecule. The energetic cost of initializing a loop is 
eo > 0, and the configuration of an open loop of size i does not further affect its energy. 
The number of states of a long loop of length £ ^ 1 is given by the number of random 
walks of 2i steps which return to their starting point: 

m - J,- (74) 

Here s is a geometrical constant which depends on the microscopic details of the 
molecule, while the universal exponent b depends only on space dimension and on the 
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existence of long range interactions in the molecule such as self-avoiding interactions: 
for a non-self- avoiding loop in d dimensions h = d/2, while self-avoiding interactions, 
both within the loop and between the loop and the rest of the molecule, were shown 



to increase the value of b to approximately 2.11 in d = 3 dimensions |48l . |49|. The 
value of this exponent has received much attention, since it determines the order of the 
transition: for 1 < b < 2 the transition is second order, while b > 2 leads to a first order 
transition. 

In recent years, with the advent of single molecule experiments, direct 



measurements of the dynamics of denatured segments became possible [50|, iMl- In 
particular, the state of a single tagged base pair can be followed using fluorescence 
correlation spectroscopy, whereby fluorescence occurs as long as the base pair is open and 
is quenched when it is closed. Such experimental developments have lead to a theoretical 
effort to study the dynamics of denaturation using the Poland-Scheraga model 
These studies consider dynamics which obey detailed balance with respect to the Poland- 
Scheraga free energy: if w±{i) are the rates with which a loop of length i changes its 
length by ±1, then w^{i)/w^{i + 1) = se~'^''{jj-r)^. At the transition temperature 



Tm, open loops are sparse. They rarely coalesce or split up since eo — lO/csTm |52 
Therefore, to a good approximation, the dynamics of a single loop may be considered 
independently of that of other loops. From these considerations one may conclude that 
at the melting temperature, the loop-length probability distribution evolves according 
to a master equation which, when the loop size is large, approaches the Fokker-Planck 
equation ([T]) (where x is the loop size i). Note that the large value of eo implies that 
once the length of a loop shrinks to zero it does not reappear in the same position for a 
long period of time. Therefore, an absorbing boundary condition at £ = is appropriate 
for the study of the dynamics of denatured loops. 

The fluorescence correlation function which can be measured in experiments is 
related to the probability that an unbound loop remains open after time t Jjol, 50 1. At 
late times, this survival probability is given by 

oo L oo 

S{t) = J2 Pi^, t) = Y. ^) + E ^) = ^1 W + (75) 
e=i e=i i=L+i 

where Si{t) and 5*2 (t) correspond to the contributions to the sum from small and 

large loops, respectively. Here, L 1 is a constant. Using the scaling form 

for the probability distribution, it is easy to evaluate the two sums and find that 

Slit) ~ t-/3-(fe+i)/2^ -vvhile S2{t) ~ t-/^-{''-i)/2^ where (5 depends on the initial condition, 

as discussed below. Therefore S2{t) dominates the sum and 

S{t) ~ (76) 

When the loop is allowed to fluctuate freely, the probability of selecting an initial 
loop of length £o in the steady state is 

Po(4)~4P*(4)~4^'"'\ (77) 
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which is normahzable in the case of DNA where b > 2 This is the natural 

experimentally-relevant initial condition when one probes the state of a base pair 
(whether it is bound or not) at random. According to the definition of the parameter a 
for the case of absorbing boundary conditions, it corresponds to a = — 1 (see Eqs. (ITSll 
and (|571) ). Thus, for the relevant initial condition one has from Eq. (fT6ll /3 = —1/2, 
yielding S{t) ~ t-{''-2)/2^ ^^^g previously obtained in 0] using a different method. 

A different possible experimental protocol is obtained when one forces one end of 
the loop to be on a particular site. In this case, there is no need for the factor of io in 
Eq. ([77D [7]. The initial condition is then Po(4) = ^*(4), yielding a = and /3 = 0. 
Therefore, the survival probability decays as S{t) ~ t-C'-i)/^^ once again in accordance 
with f^. 

Finally, the case of a localized initial condition, namely starting from a loop of a 
given length, has been considered by |8j] . This case is far harder to realize experimentally. 
In our approach, this initial condition corresponds to a = oo, which leads to /3 = 1 and 
a different behavior of the survival probability, S{t) ~ 

The conclusion from this discussion is that since the initial condition selects the 
value of the scaling exponent /3, it may affect all correlation functions which can be 
measured experimentally. Therefore, when analyzing experiments which measure the 
dynamics of denaturing DNA loops, one must carefully take into account the appropriate 
initial condition which is relevant to the experiment. 

5.4- Other systems described by Eq. ^ 

In light of its simplicity, it is not surprising that Eq. ([1]) arises in many different contexts. 
We now briefly review some of the problems described by this equation. This review, 
which is far from being exhaustive, is included to indicate the variety of problems to 
which the results obtained in this paper may apply. The physical implications of our 
results to these systems have so far not been worked out. 

(i) We have considered so far only one-dimensional problems of diffusion in a 
logarithmic potential. In fact, as long as the problem is spherically symmetric, 
diffusion in a logarithmic potential in any dimension leads to an equation of the 



form ([2])-([3]) for the diffusion in the radial direction [53|. In this case, the parameter 
b depends on the spatial dimension. A similar equation results when considering a 
spherically symmetric convection-diffusion equation in two dimensions with a sink 



or source at the origin 54 



(ii) The one- dimensional diffusion equation in an attractive logarithmic potential can 



be mapped, as we show in Appendix A to the (imaginary time) Schrodinger 



equation which describes a quantum mechanical particle in a repulsive inverse 
square potential Vs{\x\ ^ 1) ~ ■y/x^ (where the coupling constant 7 is related 
to b, see Eqs. (lA.4l) -( |Ar7l) ). The quantum inverse square potential has drawn much 



attention over the years (for references, see for example 55|, |56|). Although the 



questions we address in the present study are motivated by problems of diffusion. 
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the scaling solution we have found above is valid also in the corresponding quantum 
system. It would be interesting to understand its implications in the context of 
quantum mechanics. 

(iii) Models of gases with long-range interactions exhibit slow relaxations towards 
equilibrium. One approach to study these slow relaxations is to examine the 
evolution towards equilibrium of a single tagged particle inside an equilibrated 
gas of this type. In several models it has been established that the kinetic equation 
which describes the relaxation of the tagged-particle momentum distribution can be 
transformed to a Fokker-Planck equation with the asymptotic form ([T]), from which 



the time dependence of different correlation functions can be calculated [l2|, [iSj (for 
a review see Sec. 5.2.3] and references therein). 

(iv) An equation of type ([1]) was encountered in the dynamics of a two dimensional XY 
model below the Kosterlitz-Thouless transition. In [sS^, it has been shown that 
this equation can describe the annihilation of vortex-antivortex pairs during the 
relaxation to equilibrium after a quench from high temperatures. 

(v) In the study of Barkhausen noise, this equation is used to derive the distribution 



of magnetization jumps within the mean-field ABBM model 57|, Sect. IV, B]. 



(vi) In a biological context, a discrete-time version of Eq. ([T]) was sug gest ed as a model 



for the dynamics of sleep-wake transitions during a night's sleep [15 . 

(vii) Many studies of Eq. ([1]) were motivated not by specific physical phenomena, but 
by interesting mathematical features of the equation. These include studies of the 
persistence exponents for a diffusion described by Eq. ([1]), which are found to depend 



on the dimensionless coupling constant h 2l|, |53|, |58|] ; an examination of the effect of 



noise on evolution equations such as ([T]) which give rise to finite time singularities 
22j; and an examination of the relation between the tails of stationary distributions 
of Markov processes and power-law decay of correlations in the dynamics 59 . 

6. Conclusion 

In this paper we considered the late-time scaling behavior of a particle diffusing in 
a potential with logarithmic tails, focusing solely on the trapping case in which the 
probability distribution relaxes to a normalizable steady state. By concentrating on the 
deviation from equilibrium (i.e., the difference between the solution and the steady 



state), we have generalized the scaling solution which in 18l-l20| was obtained for 
localized initial conditions to any initial condition. 

The scaling solution to this, rather simple, linear diffusion problem contains several 
surprises. The first is that at small values of where the diffusive (x ~ -y/t) scaling 
regime is invalid, the solution is given by a different scaling function. Thus, to leading 
order in t, the full solution on the entire real axis is given by the simple scaling form ( ITOl) . 
With this new result it is easy to compute the time dependence of many correlators, 
even of functions which are concentrated around the origin (e.g., (l/x(t))). The utility 
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of this scaling form was demonstrated in the calculation of current correlations in the 
zero-range process presented in Sec. 15.21 

Another surprising aspect of the solution is that the Fokker-Planck equation (I2])-(|3]) 
has incomplete scaling solutions, i.e., solutions in which the scaling exponents cannot be 
determined from dimensional analysis. Moreover, these scaling exponents depend on the 
initial condition via a selection mechanism which is similar in many of its details to the 
marginal stability mechanism which governs selection in problems of fronts propagating 
into an unstable state. Since our system is not spatially homogeneous, the standard 
techniques which are employed in the study of the selection of propagating fronts (most 
notably Fourier analysis) are inapplicable. However, as the diffusion equation is linear, 
it can be solved exactly and the selection mechanism can be proven rigorously. We hope 
that the similarities and differences between the problem we have studied here and the 
selection in propagating fronts might shed light on the mathematical structure which 
underlies the selection mechanism. 

Beyond their intriguing and surprising mathematical properties, these scaling 
solutions have considerable utility for a large variety of physical problems which are 
mathematically equivalent to the diffusion equation in a logarithmic potential. We 
demonstrated the applicability of our results to three examples: an experimental 
protocol was suggested, in which cold atoms in an optical lattice are "quenched" from one 
value of the diffusion constant to another, which should exhibit a relaxation that depends 
on its initial steady-state; two-time current correlations in the steady-state of a system 
undergoing a non-equilibrium real-space condensation transition were calculated; and 
it was demonstrated that initial conditions are important when analyzing experimental 
data of the dynamics of denaturing DNA loops. It would be interesting to examine how 
the dependence of the scaling form on the initial condition might be manifested in other 
systems governed by the diffusion equation ([1]). A particularly intriguing question is the 
significance of such scaling solutions in the problem of a quantum mechanical particle 
in an inverse square potential. 
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Appendix A. Derivation of the scaling solution 

In this appendix we solve exactly the Fokker-Planck equation (l2])-([3l) and calculate 
its long-time asymptotic form. We proceed by performing the calculation only for 
symmetric potentials which are exactly equal to a logarithm for large enough x. The 
scaling argument of Sec. I2.6l implies that the long-time asymptotics we thus obtain hold 
for any potential with the asymptotic form ([2]). The calculation presented below is 



Diffusion in a logarithmic potential: scaling and selection 



31 



based on the methods of j3| and 58 1. The case of a = has recently been analyzed in 



a similar fashion in 



60 



Appendix A.l. Mapping to a Schrddinger equation 

Consider a particle diffusing under the influence of a symmetric potential 



V{x) 



V{x) for |x| < xq 
61og(|x|) for |x| > Xq. 



(A.l) 



for some Xq > 0, where V{x) is some symmetric potential. By a proper rescaling of x 
and t, it is always possible to set the the threshold xq = 1. We further assume that the 
potential is measured in units of temperature (i.e., ksT = 1). From now on we denote 
by a tilde any quantity in the region |x| < Xq = 1. 
The corresponding Fokker- Planck equation reads 

dP{x,t) d r , ^, ,1 d^P{x,t) 



dt 



dx 



V'{x)P{x,t) 



dx"^ 



-|e where Z = J e~ 



(A.2) 
^dx. We 



Its normalized stationary solution is given by P* {x 
wish to solve the general initial value problem defined by this equation together with an 
initial condition Po{x) = P{x, 0). As discussed in Sec. Ej by considering deviations from 
the equilibrium distribution, we may, without loss of generality, restrict our discussion 
to initial conditions with zero normalization. We therefore assume from now on that 



j Po{x)dx = 0. 



(A.3) 



To solve the initial value problem defined by such initial conditions we transform 
Eq. (]A.2|) into an imaginary-time Schrodinger equation via the transformation 23 1 



P{x,t) 



-V{x)/2 



i){x,t). 



The resulting equation for the "wavef unction" ip is 

dip{x,t) d'^ilj{x,t) 



dt dx"^ 
with the Schrodinger potential 

For the potential flA.ll) this gives 



K(x)^(x,t) 



VAx] 



with the constant 7 



Vs{x) for |x| < 1 
for \x\ > 1 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



b(b 
2 V2 



+ 1). For large x this equation describes a quantum particle 



moving in a repulsive inverse square potential. 
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Appendix A. 2. Eigenfunction representation of the solution 

By separation of variables ip{x,t) = ipk{x)T{t) we find 

Tk{t) = e-^'\ k>0 (A.8) 
which yields the time-independent Schrodinger equation 

- Vs{x)Mx) + = -k^M^)- (A.9) 

The zero energy (i.e., k = 0) eigenfunction, which corresponds to the steady-state 
solution of the Fokker-Planck equation, is 

V^*(x)=V^,=o(x) = ^e-^(^)/2 (A.IO) 

where the normalization ensures that J ip*(xydx = 1. The rest of the 

eigenfunctions can be chosen to be either even or odd, since Vs{x) is a symmetric 
potential. Denote the even eigenfunctions by tp^^ki^) and the odd by tp_^k{x), with 
A; > 0. These eigenfunctions are 

{ij±,k{x), \x\ < 1 

\A^\[c±,j{k)Jpik\x\) + c±^Y{k)Yp{k\x\)], x >1 
±y^\[c±,jik)Jp{k\x\) + c±,Y{k)Yp{k\x\)], x < -I 

(A.ll) 

where Jp and Yp are Bessel functions of the first and second kind of order p = (6 + l)/2, 
and ip±,k{x) are the even and odd eigenfunctions of the potential Vs{x). We choose to 
normalize ip±,k by demanding '?/'±,fc(l) = 1- The constants c±j and c±,y can be found 
by proper continuity requirements on the eigenfunctions at x = 1. Continuity of the 
probability P{x, t) and of the probability current J = VP + = {^"^ + |^)^~^'^^ (^^^ 
Eq. (1A.2I) ) dictate that for small e 

^±,fc(l + e)=^±,fc(l-e) + 0(e) 

^±,.(1 + ^) + l = ^±..(1 - ^) + ^ + 0{e). (A.12) 
This in turn gives for c±^j and c±^y 

c±Ak) = - v^{k)]Yp{k) + kY;{k)), 

c±Ak) = f - '-±^]Uk) - kJ'^ik)), (A.13) 

where we have defined v±{k) = -j/;^ ^(1) = v±^q + f±,2^^ + . . . (below we show that this 
series indeed contains only even powers of k). The overall normalization c±{k), chosen 
so that for large x the eigenfunctions satisfy 'ip±,kix) ~ n^^^'^ sm{kx — (l)±^k) with some 
phases 0±,a;, is 

c±ik) = f „ .J''^,^ (A.14) 



^c±Ak? + c±Ak? 
This choice of normalization guarantees the completeness relation 

/•oo 

i)*{x)i)*{x')+ I i)+^k{x)^+,k{x')+^-,k{x)^-,kW) dk = 6{x-x').{A.15) 
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Using these eigenfunctions and the definition (lA.4p . we can write down the solution 
to the original Fokker-Planck equation flA.2|) for any initial condition Pq{x). Denoting 



this solution by P(x, t|Po); "we have 

roo v(x) roo 

P(x,t|Po) = Ve-^(^)/M dka±{k)ij±,k{x)e-''"' + ^— Po{xo)dxo (A.16) 

_l_ JO ^ J~oo 

where the amplitudes a±{k) are given by the projection of the initial condition on the 
appropriate eigenfunctions 

/•oo 

a±(fc)= / rfxoPo(xo)e^(^°)/V±,fc(a;o). (A.17) 



The second term in the rhs of Eq. flA.16P is obtained by projecting Po(x) on iIj*{x), i.e., 
substituting ( lA.lOp into the expression e~^^^^^'^ip*{x) J cLxq Po(xo)e^^^°^/^'?/'*(a;o). For the 
zero-normalization initial condition ( IA.3|) which we consider, this term vanishes. 



Appendix A. 3. Eigenfunctions and amplitudes at small k 

We are interested in the long-time behavior of the solution. The e~^^* term in the first 
integral of ( 1A.16I) implies that when t ^ 1 only small values of k will contribute to 
the integral. We are therefore led to investigate the small k behavior of the amplitudes 
a±{k) (which according to the definition (1A.17I) may depend on the initial condition). 

First, let us examine the small k asymptotics of the constants c±^j{k), c±y{k) and 
c±{k). For small /c, the eigenfunctions in the region — 1 < x < 1 can be expanded as a 
power series 

i'±,k{^) = ^±fl{x) + k^h±,2ix) + k^h±^i{x) + .... (A.18) 
This expansion is uniform in x in this region, and includes only even powers of k as the 
eigenvalue problem (1A.9P is even in k. The zeroth order terms are 

tlj*{x) 



^ r{x)j;e^^y)dy 

Here, ip*{x) is given by flA.lOp . from which we can deduce that f+_o = = ~^~2^j 

and similarly t>_ q = —^-^ + ji-^^v(y)^y - substituting the expansion (]A.18|) in the 



Schrodinger equation (lA.SP and continuing the perturbative calculation to the next 



order, it can also be shown that = ^+2(1) = ~ /-i ^T"*^^' which, together with 
e~^dx = l/{b — 1), yields Z = e~^dx = 2/{b — 1) — 2t>+^2- Substituting these 
in dAH]) and (lAli]) we find 

c+,jik) = -ZT{p)Q'"\l + Oik% 

-TX {k\P^_^ ^/,2a 



-1 (k 
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and similarly c^^jik) ~ k~P, c^y{k) ~ k^ and c^{k) ~ A;''+^/^ (the coefficients of the 
latter three are omitted because they will not be used below). With these, together with 
the known asymptotics of the Bessel functions 2^, we may rewrite the eigenfunctions 
f lA.ll|) for X > 1 and A; < 1 as 



i)±,k{x > 1) 




^Ukx) + 



k\^^^ kx 




Yp{kx) 



{kx) 



r(p+i)25 



l + 0{{kx) 



kz 



X 



ZV{p)2-2- 



l + 0{{kx) 



(A.21) 



In order to study the amplitudes (]A.17|) . we must make some assumptions about 
the initial condition Pq{x). Below we assume that Po{x) is asymptotically symmetric 
for large i.e., 



Po(a; > 11 



-X < -1) ~ P*{x) ■ A\x\ 



(A.22) 



This assumption is made solely for notational simplicity. In general, one could have 



Pnix — 7- ±Oo) 



P*{x) ■ A±\x\ 



The calculation which we present below can be 
repeated for this more general case, resulting in different scaling behaviors for positive 
and negative x's, in which case only the smaller of a_|_ and a_ dominates the eventual 
long-time behavior. We further assume, without loss of generality, that for all > 1, 
not just for large x, the initial condition is already close to its asymptotic form, i.e., 
Po(kl > 1) ~ P*{x) ■ (one can rescale x and t to ensure that this is the case; 

note that such a rescaling entails a redefinition of V, ip and v±{k)). 

When the initial conditions are asymptotically symmetric, the small k behavior of 
a±{k) is determined as follows. Separating the integration in flA.171) to three integrals 
and substituting Equations ( lA.ll) and ( lA.lip . we can write 



with 



a+{k) = c+{k) y2h{k) + 2h{k) + J+,3(A;) 
a_(A;) = c_(fc)/_,3(A:), 

/oo 
dx Cj^^,j{k)J p{kx)x^ Pq{x) 

/oo 
dxc+^Y{k)Yp{kx)xPPQ{x) 

I±M= f' dx^±4x)e^^^y'Po{x) 



(A.23) 



(A.24) 



(note that Po{x) need not be symmetric for —1 < a; < 1). Changing the integration 
variable in the first integral to z = kx and substituting Equations ( 1A.20I) . ( 1A.22I) and 

P*(x) 



^e-^(^), yields 



h{k) ^ -Ar(p)2^-2 . r ■ 



(A.25) 
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The latter integral converges to a finite value when — )■ if (1 — 6 — 2a)/2 + p > —1, 
or equivalently if a < 2; otherwise it diverges with k. Its asymptotic behavior is given, 
to leading order in A;, by 

' 2i-P-T(l-§) 



r(^i+§±i) 

2-(b+l)/2 



when a < 2 



P^fc+S) log k when a = 2 . (A. 26) 



1 U2-a 



when a > 2 



(a-2)2Pr(p+l) 

To evaluate hik), for all fcx ^ 1 we can approximate c+^Y{k)Yp{kx) = x~P[l + 0{{kx)'^)] 
(using Eq. (IA.20p and the known asymptotics of the Bessel function). Fixing an e ^ 1, 
the integral is evaluated as 

pe/k /"OO 

l2(k) = / Poix)dx + / c+y{k)Yp{kx)x''Po{x)dx + 0{k'^) = 

Jl Je/k 

/•OO /"OO 

= / Po{x)dx+ [c+y {k)Yp{kx)xP -l]Po{x)dx + 0{k'^). {A.27) 

Jl Je/k 

Changing once again the integration variable to z = kx reveals that 

/OO 
Po{x)dx + O{k^,k''-^''-^). (A.28) 

The integrals I±,3{k) are evaluated using the expansion ( ]A.18p - fA.19|) : 
/+,3= / Po{x)dx + O{e) 



^-'3= / ri ^, ,' Poix)dx + Oik'). (A.29) 

Combining (1X251) - (IA:26D . ^K28\f and (1X29]) into Equations flA^Oj) and (1X231) . and 
remembering that Po(x)dx = 0, finally yields, to leading order, 

a±{k) ~ C±r± when a ^ 2, (A.30) 

with 

^1 J f ^ + l-(2-a) if a<2 

._ = - + !, and -. = ||^i ^2 " ^^'^^^ 

The constants C-t are non-universal (i.e., they depend on the full forms of the potential 
V{x) and of the initial condition Pq{x)), except when a < 2, where 

C.ia < 2) = I . ^^^^ (A.32) 

When a = 2, a_{k) is still given by Eqs. (]A.30p -( fA.3ip . but the expression for a+{k) is 
replaced with 

A 2"^/^ b 

a+{k)^--- ■ k-2+' log k when a = 2. (A.33) 
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Appendix A. 4- Late-time scaling solutions 

Once the asymptotic forms of ip±,k{x) and a±{k) for small k (Eqs. ( ]A.19|) . ( 1A.21I) . 
( ]A.30|) -( fA.33l) ) are known, they can be substituted into equation ( 1A.16I) . When a 7^ 2, 
changing the integration variable to g = kt^^'^ yields in the region |x| < 1 

POO 

dka+{k)ij+4x)e-^'' = t'^'^ / dqa+{qt-^'^)c+{qt-^/^)i)^^^-^/2{x)e-'i' = 



6+2^ I 

c+t — —ij*{x) 



(A.34) 



A similar calculation for the odd eigenfunctions shows that their contribution is 
negligible in comparison with (lA.34p for all |x| < 1. In the region |x| > 1, we similarly 



have 

dk a±{k)i/j±^k{x)e 



c+r 



4X2 



r(^)2t+i 



X 



X 



t 



2(6+1) 



X) 

'dq- ^'"^ ^' / gt+-±-ie-^'dgJ (^l + 0(xV^)j. 

(A.35) 

When a = 2, equations (1A.34I) and (]A.35P have a similar form but are multiplied by an 
overall logt correction factor. 

As long as \x\ <^ t^/^, the higher order terms in Eqs. flA.341) and flA.35P can be 
dropped. Using the identity q'^e~'^^dq = r(^^)/2 then leads to 



P{x < Vt,t\Po 
P*{x) 



b+a-l 
t 2 



4(fe+l) _^ 



Z (b+a-l) 



C ( 



b+1 

t 2 log t 



Z t 



t 2 



where 



C 



Ar(i-f) 

A 



Z t 



when a < 2 
when a = 2 



when a < 2 
when a = 2 
when a > 2 



(A.36) 



(A.37) 



and C is non-universal when a > 2 (compare with Eqs. (ITU]) -( IT^ and (IT^ - (IT5]) ). Note 
that Eq. flA.361) holds even in the region |x| < 1, where the potential is not logarithmic. 

When X = 0(t^/^), the higher order terms in Eq. flA.35p cannot be neglected. They 
are taken into account by leaving the Bessel functions in Eq. (]A.21|) unexpanded when 
substituting in Eq. f lA.16|) . Changing the integration variable once again to g = kt^^'^, 
and substituting u = xt'^^"^, yields when a ^ 2 



dk ilj±^k{x)a±{k)e' 



C+vflH- 
V2 



v± + l 



Jp{qu) + 



6-1 



Zr2(p)2^- 



r^Yp{qu)t- 



b-l 
2 



dq, 
(A.38) 
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and a similar expression with a logt correction when a = 2. The second term in the 
square brackets is neghgible at late times. Using the identity 0] 



q^Jp{qu)e dq 



r(i + p) 



1 + p + /i 



p 



(A.39) 



together with 

P(Mti/^t|Po 



P*{uty^ 



we finally arrive at the scaling solution 



,6+1 rp / 6+0+1. 6+3. 1 ^-^ 



^6+lg-2i- . ^■ 



-1 



when a < 2 
when a = 2 
when a > 2 



(A.40) 



where the constant C is the same as in (1A.37P (compare with Eqs. (ITU]) - (IT^ and ( IT^ - 

(USD)- 

We have repeated the calculation of this appendix also for the case considered in 
Sec. Hlof absorbing boundary conditions at the origin. This lengthy but straightforward 
calculation, which we do not present here, recovers Eqs. ( 158|) -( 16T|) and yields the 
logarithmic corrections when a = 2. In particular, the calculation reveals that Eq. 
(IA.40p holds regardless of the boundary condition. 



Appendix B. Calculation of the sum in Eq. (I70p 



In this appendix, we calculate the current correlation function discussed in Sec. 15. 2 1 using 
the scaling solution (fTOj) . We begin by splitting the sum in Eq. (170|) into three terms. 



C{t) = Si{t) + S2{t) + S3{t), 

where we define 

Slit) 

S2{t) 



(B.l) 



[Pim,t\Po) - P*{m) 

m=ni (t)+l 
ni{t) 

^P{m,t\Po) - P*{m) 



(w(m] 



w[m) 



(B.2) 



m=l 

oo 



P*(m) 



ni{t) 

n 



P(mA\Pn) - P*(m) 



m=0 



m=ni {t)+l 

Here, ni(t) is chosen to satisfy <^ ni{t) <^ t^^"^, and in the last equality we have 

used the normalization condition '^mP*ijn) = J2mP(™^t\Po) = 1. 

As discussed in Sec. 15.21 at large times, the terms in the square brackets in (1B.2P 
can be replaced by the scaling solution ffTOj) and f|T6|) with a = 1. In the first sum, the 
square brackets are replaced with the large-x scaling function, yielding 



Slit) 



E 



P*iuVi)Ct-'^'fy2iu)] (wiuVi) - 1 
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Ch 6+1 







°° ^ /1 + 6 + a 6 + 3 u''\ ^ 

i-Ti ; ; du 

\ 2 ' 2 ' 4 



7rr(i±^) .+1 

Z2^r2(|) 

(B.3) 



where we have substituted the asymptotic form of P*{n) (168!) and the value of the 
constant C which is given in (fTSl) . 

A similar calculation is carried out for 5*2 and 5*3, this time using the small-x scaling 
function. We now show that although 82,83 ~ t'^^"^, the two sums cancel each other 
to leading order in t. To see this, substitute (ITUl) in flB.2p . and use fIBSl) to deduce that 
P*{n)w{n) = P*{n — 1). Combining these gives 



82{t) ~ J2t-''^^P*{ 



■mjwim) 



m=l 
ni(t)-l 



C3 + C 



m 



J2 t-'/'p* 



m=0 
ni{t) 



C. + C 



Cs + C 



'm + 1 



t 



(B.4) 



m 



b+l 



where C3 is given in (jH]). Therefore, 

ni — 1 

82{t) + 53(t) ^ 5^ Crt-^P*(m)[(m + 1)'+^ - m'+^] - 



m=0 



t-^P*{n,{t)) 



C, + C- 



t 



(B.5) 



6+1 



from which it can be shown that, to this order, \82{t) + 8s{t)\ ^ t~~. Using a 
perturbative expansion similar to ( ]A.18|1 . it can be shown that the contribution from 
higher order corrections to the scaling form are also negligible. 

Adding the three contributions together, we find that the asymptotic decay of the 
correlation function is 



C{t)^8^{t) 



,_6+l 
■t 2 , 



(B.6) 
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